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ABSTRACT 

We present 1.4 and 5 GHz observations taken with the Very Large Array (VLA), and 
observations at 250 GHz obtained with the Max-Planck millimeter bolometer (MAMBO) at the 
IRAM 30 m telescope, of ten optically selected Quasi-stellar Objects (QSOs) at 5.0 < z < 6.28. 
Four sources are detected at 1.4 GHz. Two of the sources have rest- frame 1.4 GHz luminosity 
densities > 5.0 x 10 26 W Hz -1 , placing them in the regime of radio-loud QSOs. Both of these 
sources are also detected at 5 GHz. These results are roughly consistent with there being no 
evolution of the radio-loud QSO fraction out to z ~ 6. 

Three sources have been detected at 250 GHz or 350 GHz by these, and previous, 
observations. The (sub)mm flux densities for these three sources are much larger than their 
1.4 GHz flux densities. The rapidly rising spectra into the (rest frame) far IR argue that 
the observed mm emission is likely thermal emission from warm dust, although more exotic 
possibilities cannot be precluded. The implied IR luminosities are between 10 12 and 10 13 L©. 
For J030 1+0020 the radio continuum emission is clearly above that expected for a star forming 
galaxy based on the radio-FIR correlation. In this case it seems likely that the radio emission 
relates to the AGN. For J0756+4104 the radio emission is within the range expected for a star 
forming galaxy, while for J1044-0125 the radio upper limit is at least consistent with a star 
forming galaxy. If the dust is heated by star formation, the implied massive star formation rates 
are between 200 and 1000 M Q year" 1 . 

We do not detect radio emission from the reported X-ray jet associated with J1306+0356. 
The lack of radio emission implies that the magnetic field is well below typical equipartition 
values in powerful radio jets, or that particle acceleration ceased between 10 6 and 10 7 years ago, 
or that the X-ray emission is not inverse Compton emission from a jet related to J1306+0356. 

The highest redshift source in our sample (J1030+0524 at z = 6.28) is not detected at 1.4 
or 250 GHz, but four fairly bright radio sources (Si.4gh z > 0.2mJy) are detected in a 2' field 
centered on the QSO, including an edge-brightened (TRIP) double radio source with an extent 
of about 1'. A similar over-density of radio sources is seen in the field of the highest redshift 
QSO J1148+5251. We speculate that these over-densities of radio sources may indicate clusters 
along the lines-of-sight, in which case gravitational lensing by the cluster could magnify the 
QSO emission by a factor 2 or so without giving rise to arcsecond-scale distortions in the optical 
images of the QSOs. 
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Subject headings: dust: QSOs, galaxies — radio continuum: QSOs, galaxies — infrared: QSOs, 
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1. Introduction 

Surveys such as the Sloan Digital Sky Survey (SDSS) (York et al. 2000) have generated large samples 
of QSOs out to z^6. Optical spectroscopy of two of the highest z QSOs revealed, for the first time, 
evidence for broad, optically thick absorption regions in the (rest frame) far UV, as expected for the neutral 
intergalactic medium (Gunn & Peterson 1965; Becker et al. 2001; Djorgovski et al. 2001; Pentericci et al. 
2002; Fan et al. 2003), suggesting that we are probing into the epoch of reionization, the edge of the "dark 
ages" when the first stars and massive black holes were formed. Assuming that the intrinsic continuum 
spectra of the z > 5 QSOs are similar to those of nearby QSOs as compiled by Elvis et al. (1994), and that 
these QSOs are not lensed, Eddington limit arguments result in lower limits of several times 10 9 M© for the 
masses of the black holes in these high redshift objects. 

Study of the dynamics of stars and gas in the nuclear regions of nearby galaxies has led to two 
remarkable discoveries: (i) the overwhelming majority of spheroidal galaxies in the nearby universe contain 
massive black holes and (ii) there is a correlation between the black hole mass and the velocity dispersion of 
the stars in the spheroid (Magorrian et al. 1998, Ferrarese & Merritt 2000, Gebhardt et al. 2000, Trcmainc 
et al. 2002). Assuming that these relations also apply in the high redshift universe (Shields et al. 2002), 
implies that the highest redshift QSOs are associated with massive galaxies (> 10 11 M ). Furthermore, 
strong metal emission lines in these QSOs suggests that their environments have been quickly enriched, 
possibly through starburst activity (Fan et al. 2001). 

These fundamental results allow us to investigate the nature of these earliest objects, and in particular 
to probe the relationship between massive black hole and spheroidal galaxy formation. To address this and 
other related questions, and in general, to study the cosmic evolution of the radio-to-optical spectra of 
QSOs, we have undertaken an extensive observational program from radio through millimeter wavelengths 
of high redshift QSOs which includes: searches for emission from warm dust at millimeter wavelengths from 
a large sample of z > 1.8 QSOs (Carilli et al. 2001a, Omont et al. 2001, 2003), high resolution imaging 
at cm wavelengths of the non-thermal radio continuum emission from these sources (Carilli et al. 2001a, 
b), and observations of the CO line emission from selected sources with large infrared luminosities (Carilli, 
Menten & Yun 1999; Carilli et al. 1999, 2002a, b; Cox et al. 2002; Beelen et al. in prep.). 

Millimeter-continuum observations of high redshift QSO's show that about 30% of the sources are 
detected in surveys with flux density limits of 1 to 2 mJy at 250 GHz (Omont et al. 2001; Carilli et al. 
2001b; Isaak et al. 2002). Multifrcquency studies of the rest-frame radio through far IR (FIR) spectral 
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energy distributions (SEDs) of these sources reveal that the mm emission is thermal emission from warm 
dust (Bcnford et al. 1999), and in many of the sources the cm-to-mm flux density ratios are consistent with 
the radio-to-FIR correlation for star forming galaxies (Carilli et al. 2001b). Furthermore, searches for CO 
line emission from FIR luminous QSOs have resulted in the detection of molecular line emission indicating 
the presence of large gas reservoirs (~ 10 11 M Q ; Cox et al. 2002). These detections lead some authors 
to conclude that active star formation is inevitable (Omont et al. 2001) in these FIR luminous, CO-rich 
quasars. Andreani et al. (2002) reached a similar conclusion for lower redshift QSOs, although for both low 
and high redshift samples the case for star formation is by no means proven. 

In this paper we extend our radio study of distant QSOs to the highest redshifts (z > 5). We present 
observations at 1.4, 5, and 250 GHz. These observations are an order of magnitude, or more, more sensitive 
than survey observations, such as FIRST and NVSS (Becker, White, & Helfand 1995; Condon et al. 1998). 
We assume a concordance cosmology (fl m = 0.3, il\ = 0.7) and H = 65 km sec -1 Mpc -1 throughout. 

2. Observations 

We observed with the VLA all z > 5 QSOs known at the time of observation, except J2245+0024 
(Sharp, McMahon, Irwin & Hodgkin 2001), which was outside our LST range. The sources are listed in 
Table 1. Most of the sources come from the SDSS, and hence are luminous QSOs (Mb < —26), with the 
exception of J0301+0020, which is a low luminosity QSO discovered by Stern et al. (2000b). 

VLA observations at 1.4 GHz were made in the A configuration (max. baseline = 30 km), with a total 
bandwidth of 100 MHz and two polarizations. Each source was observed for about two hours. In addition, 
J0756+4104 was observed for 2 hours at 1.4 GHz in B configuration (max. baseline = 10 km). Standard 
phase and amplitude calibration was applied, and all sources were self-calibrated using field sources. The 
absolute flux density scale was set with observations of either 3C48 or 3C286. 

The final images were generated using the wide field imaging (Cotton 1999; Bridle & Schwab 1999) and 
deconvolution capabilities of the AIPS task IMAGR. The theoretical rms noise (a) value corresponding to 
2 hours of observing in continuum mode at 1.4 GHz is 16 /iJy, and in most of the maps presented here this 
sensitivity is roughly achieved. The noise level in one source (J1204-0021) was higher (35 /xJy), possibly 
due to some low level terrestrial interference. We also include a 2% uncertainty in absolute flux density 
scale determination. The Gaussian restoring CLEAN beam Full Width at Half Maximum (FWHM) was 
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typically ~ l'/5 for the A configuration observations. 

Two of the sources, J0836+0054 and J0913+5919, were observed at 5 GHz on August 17, 2002 in B 
configuration for about 20 minutes, achieving an rms sensitivity of about 60/iJy 1 . The Gaussian restoring 
CLEAN beam (FWHM) was ~ 1'.'5 , matching the resolution of the 1.4GHz observations. 

Observations at 250 GHz were made using MAMBO (Kreysa et al. 1999) at the IRAM 30m telescope 
on Pico Vclcta in Spain, during the winter of 2001-2002 within dynamically scheduled, pooled observations, 
except for J0301+0020, which was reported on by Bertoldi & Cox (2002). MAMBO is a 37 element 
bolometer array with an effective central frequency of 250 GHz for thermal sources. The beam for the feed 
horn of each bolometer is matched to the telescope beam of 10'.' 6. Observations were made in standard 
on-off mode, with 2 Hz chopping of the secondary by 50" in azimuth. The data were reduced using the 
MOPSI software package (Zylka 1998). Pointing was monitored every hour, and was found to be repeatable 
to within 2". The sky opacity was monitored every hour. Zenith optical depths were lower than 0.3 at all 
times. Gain calibration was performed using observations of Mars, Uranus, and Ceres. We estimate a 20% 
uncertainty in absolute flux density calibration based on these observations. 

The target sources were centered on the central bolometer in the array, and the temporally correlated 
variations of the sky signal (sky noise) detected in the surrounding six bolometers were subtracted from the 
central bolometer signal. The total sky + source observing time for each source was about 1 hour, leading 
to rms sensitivities of 0.5 to 1.0 mJy, depending on the weather. 

3. Results 

The results from the 1.4 and 250 GHz observations are listed in Table 1. The abbreviated source 
names (column 1), the redshifts (column 2), the absolute blue magnitude (Mb) obtained by assuming a 
flat A dominated cosmology (column 3), and the optical position (columns 4 and 5) are compiled from the 
optical discovery papers. Columns 6 and 7 give the 1.4 and 250 GHz flux densities with la error-bars, and 
3cr upper limits are given for non-detections. 

The 1.4 GHz images of these sources are shown in Figure 1. The positional uncertainty for the radio 

1 D. Rusin and collaborators also observed J0836+0054 on May 5 at 5 GHz. Their flux densities are equal 
to ours within the errors. 
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observations is given by: ag ~ ^nr 1 (Fomalont 1999), where FWHM corresponds to that of the Gaussian 
restoring beam, and SNR=Signal-to-Noise ratio of the detection. For a 3a detection this corresponds to 
0'.'5 for most of our sources. To this must be added the typical astrometric uncertainty 0>!l (Pier et al. 
2003) of the optical data, and the uncertainty in the relationship between radio and optical reference frame 
which is about 0725 (Deutsch 1999). Thus, we only consider emission within 0'.'6 of the optical position to 
be associated with the QSO. 

Fomalont et al. (2003) show that the sub-mJy source counts follow the relation: 
N(> Si. 4 ) = 0.026 x Sri' 1 arcmin^ 2 , with 1.4 GHz flux density Si. 4, in mJy. Hence, within 
0'.'6 of a given source we expect 1.5 x 10~ 4 sources with Si. 4 > 70/zJy by chance. Blain et al. (2002) show 
that there are about 2000 sources deg -2 with S250 > 1 mJy. At this flux density level we then expect 0.01 
sources by chance within the beam of the 30m telescope (ie. within 5" of the target source). 

Notes on individual sources 

J0231 0728 (z=5.41) There is a 3a unresolved source at 1.4 GHz situated l'.'O away from the 
optical position of the QSO. This is farther than the 0'.'6 positional accuracy of our measurements so 
J0231-0728 is considered a non-detection at 1.4GHz. This source is also not detected at 250 GHz. 

J0301+0020 (z=5.50) A radio continuum source with S1.4 = 73 ± 17^Jy is detected within 0'.'5 of 
the optical QSO position. The radio source is not resolved, and Gaussian fitting sets an upper limit to 
its size at 1.4 GHz of l'/6 . Bertoldi & Cox (2002) detect this QSO at 250 GHz with a flux density of 
0.87 ± 0.20 mJy. 

J0756+4104 (z=5.09) The combined A and B configuration observations show a 1.4 GHz source 
with Si. 4 = 65 ± 17^iJy within 0'.'2 of the optical QSO position. Gaussian fitting to the radio emission sets 
an upper limit to its size of 2"3 . This QSO is also detected with MAMBO with S250 = 5.5 ± 0.5 mJy. 

J0836+0054 (z=5.82) This source is clearly detected at 1.4 GHz with a flux density of 1.75 ± 0.04 
mJy. Gaussian fitting sets an upper limit to the source size of 0'.'65 at 1.4 GHz. J0836+0054 is also detected 
at 5 GHz with a flux density of S5 = 580 ± 57yuJy. The implied radio spectral index is -0.8. A second radio 
source with Si. 4 = 0.44mJy is detected 10" south of the QSO. Deep optical and near-IR imaging of this 
field suggests that this second source is associated with a lower redshift galaxy, and unrelated to the QSO 
(Rusin et al. in prep). J0836+0054 is not detected at 250 GHz with a 3a upper limit of 2.9 mJy. 

J0836+0054 was also detected in the FIRST radio survey (Becker, White & Hclfand 1995) at 1.4 GHz 
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with a flux density of 1.1±0.15 mJy. The difference in flux density between the FIRST measurement and 
our more recent measurement is significant at the 4cr level, such that the radio source appears to be variable 
on yearly time-scales. 

J0913+5919 (z=5.11) This source is clearly detected at 1.4 GHz with S1.4 = 18.95 ± 0.4 mJy. 
Gaussian fitting to the radio emission sets an upper limit to the source size of 0'.'2 at 1.4 GHz. J0913+5919 
is also clearly detected at 5 GHz with S5 = 8.1 ± 0.2 mJy, implying a spectral index for this source of -0.7. 
This source is not detected at 250 GHz with a 3a upper limit of 2.8 mJy. 

J0913+5919 was also detected in the NVSS radio survey (Condon et al. 1998) with a total flux density 
of 18.5±0.5 mJy, implying that the source is not variable on yearly time-scales. 

J1030+0524 (z=6.28) This source was the highest redshift QSO published at the time of observation. 
The source is not detected at 1.4 or 250 GHz. However, four fairly bright sources (Si. 4 > 200^iJy) are 
detected at 1.4 GHz in a 2' field centered on the QSO (Figure 2), including an edge-brightened double radio 
source ('Fanaroff- Riley Class IF; Fanaroff & Riley 1974) with an extent of about 1' and a total flux density 
of 29 mJy. 

J1044-0125 (z=5.73) This source is not detected at 1.4 GHz with a 3<r upper limit of Si. 4 < 79/zJy. 
We did not observe this source with MAMBO. Iwata et al. (2001) report a flux density at 350 GHz for 
J1044-0125 of 6.2 ± 2.0 mJy. Assuming a thermal IR spectrum typical for Ultraluminous IR galaxies leads 
to an expected 250 GHz flux density of 3.4 ± 1.1 mJy for a source at z = 5.73. Near-IR spectroscopic 
observations of this source reveal a prominent CIV absorption feature whose shape suggests that this is a 
Broad Absorption Line QSO (Maiolino et al. 2001; Goodrich et al. 2001). 

4. Analysis 

4.1. Fraction of radio loud objects 

Questions regarding both the bi-modality of the radio luminosity function of QSOs, and its evolution, 
have been investigated by numerous workers (e.g. Peacock, Miller & Longair 1986; Miller, Peacock, & Mead 
1990; Schmidt, van Gorkom, Schneider, & Gunn 1995; Stern et al. 2000; Lacy et. al 2001; Ivezic et. al. 
2002). Two definitions are generally used to demarcate radio-quiet and radio-loud QSOs. One criterion is 
based on the radio-optical ratio, R ro , of the specific fluxes at rest-frame 5GHz and 4400A (Kcllerman et 
al. 1989) where typical radio-loud sources have R ro in the range 10-1000, while radio quiet sources have 
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R ro < 1. Peacock, Miller & Longair (1986) point out that R ro can be used as a discriminating parameter 
only if the radio and optical luminosities are linearly correlated, which does not seem to be the case (Stocke 
et al. 1992). The second definition divides the sources at the 1.4 GHz luminosity density of 3 x 10 25 W 
Hz -1 (Gregg et al. 1996). Ivezic et al. (2002) find that for optically selected QSOs the two definitions are 
consistent as a consequence of selection effects in flux limited samples, so for the remainder of this paper 
we will use the Gregg et al. (1996) definition to classify a source as radio loud. 

Having the radio spectral indices of the two bright radio sources, J0836+0054 and J0913+5919 (—0.8 
and —0.7, respectively), allows us to estimate their luminosity densities at a rest frame frequency of 
1.4 GHz. For J0836+0054 the value is 5.0 x 10 26 W Hz" 1 , while that for J0913+5919 is 4.9 x 10 27 W 
Hz -1 . Both these sources are radio-loud by any definition. 

The radio spectral indices for J0756+4104 and J0301+0020 are unknown, so we calculate luminosity 
densities using the spectral index of —0.5 adopted in previous studies (eg. Stern et al. 2000, Ivezic et al. 
2002). The 1.4 GHz rest frame luminosities for these two sources are 1.2 x 10 25 and 1.6 x 10 25 W Hz -1 
respectively, placing them below the radio-loud demarcation. 

4.2. The radio-FIR correlation 

Three sources (J0756+4104, J0301+0020 and J1044-0125) are detected at (sub)mm wavelengths at 
levels at least ten times above their 1.4 GHz flux densities, implying rapidly rising spectra from cm-to-mm 
wavelengths. Note that an observing frequency of 250 GHz corresponds to a rest frequency of 1600 GHz 
(= 188 /jm) for a source at z = 5.5, such that the MAMBO observations sample the rest frame far IR 
part of the spectra. Benford et al. (1999) have performed multiwavelength observations of a number of 
mm-loud high redshift QSOs, and in every case they find that the rest frame far IR SEDs are consistent 
with a gray body spectrum characteristic of thermal emission from warm dust. In the analysis below we 
will assume that the rapidly rising cm-to-mm spectra of the three mm-loud sources in our sample imply 
thermal emission from warm dust. However, we cannot preclude more exotic explanations for the rapidly 
rising spectra, such as synchrotron self absorption at rest frame far IR wavelengths. 

One possible method for addressing the question of dust heating is using the fact that star-forming 
galaxies at low redshift follow a very tight linear relation between radio continuum and FIR luminosity 
(Condon & Yin 1990; Condon 1992; Crawford, Marr, Partridge & Strauss 1996; Miller & Owen 2001; Yun, 



-10- 



Rcddy, & Condon 2001). This correlation holds over four orders of magnitude in luminosity with only 
a factor of two scatter around linearity for galaxy samples selected in the optical, IR, and radio. While 
a general correlation is expected since the synchrotron radiation at centimeter wavelengths, and thermal 
dust emission at IR wavelengths, both relate to massive star formation, the tightness and linearity of 
the correlation remain puzzling. A further uncertainty in using the FIR-radio correlation, to assess the 
significance of star-formation in heating the dust, is that lower luminosity radio quiet QSOs at lower redshift 
also follow the standard radio-FIR correlation for star forming galaxies (Sopp & Alexander 1991). Whether 
the sources in the Sopp & Alexander sample also host active star formation remains unknown. Also, there 
is the question of whether the radio-FIR correlation holds in star forming galaxies at high redshift (Carilli 
& Yun 1999; 2000; Yun & Carilli 2002; Garrett 2001). Given these uncertainties, we feel that comparing 
centimeter and millimeter continuum luminosities is simply a consistency check for star formation, but not 
conclusive evidence for it. 

Figure 3 shows the relationship between redshift and the 250-to-1.4 GHz spectral index (cuf 5 ^) for a 
star-forming galaxy taken from the study of Carilli & Yun (2000). This model consists of the mean SED 
(plus scatter) of 17 low-redshift star-forming galaxies. The far IR part of this model SED corresponds 
roughly to a modified blackbody spectrum with a dust temperature of 50K and emissivity index of 1.5, and 
the centimeter part of the spectrum is constrained to follow the radio-to-FIR correlation for star forming 
galaxies. A point located below the curve (plus scatter) on this diagram indicates a source that is radio 
loud relative to a star-forming galaxy. In such an object it is likely that the radio emission is due to the 
AGN and not star formation (Yun et al. 2001). 

For J0301+0020 we find af 5 ® = 0.5 ± 0.1, well below the region defined by star forming galaxies. For 
J0756+4104 we find af 5 ® — 0.86 ± 0.08, which is consistent (within the errors) with the low end of the 
range for star forming galaxies. For J1044-0125 we extrapolate to 250 GHz from the measured 350 GHz 
flux density (section 3). This source was not detected at 1.4 GHz, implying af 5 ° > 0.7, which is at least 
consistent with star formation. A factor two deeper radio observations are required to test whether this 
source follows the radio-FIR relationship defined by low redshift galaxies. 

4.3. Inverse Compton emission from J1306+0356? 

X-ray observations of the source J1306+0356 (Brandt et al. 2002) suggest a possible jet-like feature 
23" from the QSO position (Schwartz 2002). Schwartz (2002) argues that the emission may be due to 
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inverse Compton (IC) scattering of the cosmic background by relativistic electrons in a jet emanating from 
the QSO. These same relativistic electrons would emit radio synchrotron radiation in the presence of a 
magnetic field, hence we have searched for radio emission from the location of the jet-like X-ray feature. 

We have not detected radio emission from the possible X-ray jet in J1306+0356 to a 3<j limit of 150 
yitJy at 1.4 GHz (after convolving to the 5" x 2" resolution corresponding to the box containing the X-ray 
'jet'). The flux density of the X-ray feature at 1 keV (2.4 x 10 17 Hz) is 8.3 x 10~ 4 yuJy, assuming a spectral 
index of —1. 

A long standing and well documented technique for deriving magnetic fields in extragalactic radio 
sources is by comparing the radio synchrotron and X-ray IC flux densities (Harris & Grindlay 1979). The 
constraint on the magnetic field strength comes from the fact that the IC X-ray emissivity is a function of 
the relativistic electron density and the energy density in the dominant ambient photon field (presumably 
the microwave background), and the synchrotron radio emissivity is a function of the relativistic electron 
density and the magnetic energy density. A recent simple parameterization of this calculation can be found 
in equation 4 in Carilli & Taylor (2002). In the case of the jet-like X-ray feature in J1306+0356, only a 
radio upper limit is available, such that we can only derive an upper limit to the magnetic field strength. 
Using equation 4 from Carilli & Taylor (2002), and the X-ray flux density and radio upper limits given 
above, we derive an upper limit of 3^G to the magnetic field in the jet-like feature in J1306+0356. Note 
that this limit assumes the X-ray emission is IC, and that the relativistic electrons have a power-law energy 
distribution of index -3 over a wide range in energy (see Section 5.3 for more details). 

5. Discussion 

5.1. Radio loud QSOs at high redshift 

There has been considerable debate in the literature about the redshift evolution of the radio-loud 
fraction of QSOs (e.g. Visnovsky et al. 1992; Schneider, van Gorkom, Schmidt & Gunn 1992; Schmidt, van 
Gorkom, Schneider, & Gunn 1995; Hooper et al. 1995; Goldschmidt et al. 1999; Stern et al. 2000). Most 
recently Ivezic et al. (2002) have investigated this question using large samples of QSOs out to z ~ 2.2 
from the SDSS. They conclude that about 8% of QSOs are radio-loud, independent of redshift. Two out of 
the ten QSOs at z > 5 observed in our study are radio-loud, which is roughly consistent (given the small 
number statistics) with the fraction seen at lower redshift, eg. the Poisson probability of seeing 2 radio-loud 
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quasars when the expectation value is 0.8 is 14%. As emphasized by Ivezic et al (2002), larger QSO samples 
at high redshift are required to separate effects related to rcdshift-luminosity biases in flux limited samples. 

The radio loud sources J0836+0054 and J0913+5919 are compact (< 1" ) and have steep spectra (-0.8 
and -0.7), classifying them as Compact Steep Spectrum (CSS) objects. CSS objects are a mixed-bag of 
source types, ranging from steep spectrum core-jets (size < few pc), to small (size ~ 1 kpc), double- lobcd 
radio galaxies, ie. compact symmetric objects (CSOs; Tzioumis et al. 2002). The CSOs are particularly 
intriguing since they are thought to be very young radio sources, and their ages can be measured through 
VLBI observations of the proper motions of the radio hot spots, with typical measured ages between 10 3 
and 10 4 years (Polatidis & Conway 2002). Recent VLBI imaging of J0836+0054 by Frey et al (2003) at 
1.6 GHz shows a marginally resolved source at 10 mas resolution, suggesting that this source is a steep 
spectrum core-jet, and not a CSO. Similar VLBI imaging has not been performed for J0913+5919, but the 
lack of variability of this source (section 3) makes it a prime candidate for a high redshift CSO. 

In either case, the fact that the radio-loud fraction of QSOs appears to be relatively constant out to 
very high redshift, and that the two radio-loud QSOs at z > 5 discovered thus far have steep radio spectra, 
is encouraging from the perspective of studying the neutral intergalactic medium (IGM) during the epoch 
of reionization through observation of HI 21cm absorption (the '21cm Forest'). Calculations by Carilli, 
Gnedin, & Owen (2002), and Furlanetto & Loeb (2002) show that the next generation low frequency, large 
area radio telescopes, such as the LOFAR and the Square Kilometer Array, will be able to study the neutral 
IGM beyond the epoch of reionization (EOR) in H I 21cm absorption toward discrete radio sources as faint 
as a few mJy at frequencies below 200 MHz. Both J0836+0054 and J0913+5919 would be easily adequate 
for such studies if they were placed beyond the EOR. 

5.2. Thermal emission from warm dust and star formation 

Three of the ten sources in our study (J0301+0020, J0756+4104, and J1044-0125) have been detected 
at (sub)mm wavelengths with flux densities much larger than their 1.4 GHz flux densities. Again, such 
steeply rising spectra provide strong evidence that the (sub)mm emission is thermal emission from warm 
dust. 

We have considered these sources in the context of the radio-to-FIR correlation for star forming 
galaxies, as quantified in the 1.4-to-250 GHz spectral index, a™ (Figure 3). For J0301+0020 the a\ b 2 value 
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is clearly below that expected for a star forming galaxy. In this case it seems likely that the radio emission 
relates to the AGN activity. For J0756+4104 the a 25 ^ value is within the range defined by star forming 
galaxies, while for J1044-0125 the lower limit to a 25 ^ is at least consistent with a star forming galaxy. 

Given the interest in co-eval massive black hole-spheroidal galaxy formation, it is instructive to 
consider what the properties of these sources would be if the 250 GHz emission were a result of dust heated 
by a starburst. Using the relations between 250 GHz flux density, FIR luminosity, and star formation rates, 
in Omont et al. (2003) based on typical spectra of low redshift ultra-luminous infrared galaxies (i.e. Arp 
220) leads to: L Fm = 2.8 x 10 12 L Q for J0301+0020, 1.4 x 10 13 L for J0756+4104, and 7 x 10 12 L for 
J1044-0125 2 , and massive (> 5 M ) star formation rates of 200, 1000, and 500 M Q year -1 , respectively. At 
these extreme star-formation rates most of the stars in a large spheroidal galaxy could form in a dynamical 
time of 10 8 years. However, it should again be stressed that these rates assume the dust is heated by star 
formation, as opposed to being heated by the AGN itself. 

5.3. X-ray-loud, radio-quiet jets at high redshift? 

We have not detected radio emission from the reported X-ray 'jet' in J1306+0356 (Schwartz 2002). 
The implied upper limit to the magnetic field (3/iG) is more than an order of magnitude below typical 
magnetic field values in powerful radio jets. Perhaps the simplest conclusion from these results is that the 
X-ray feature is not IC emission from a jet emanating from J1306+0356. 

It is possible that such a jet could be Inverse Compton X-ray-loud and still be radio-quiet, even for 
strong magnetic fields, due to the different radiative lifetimes of the particles involved. The spectral peak 
of the CMB behaves as 1.6 x 10 11 (l+z) Hz, such that observations at 1 keV are sensitive to electrons 
with Lorentz factors, -f e ~ 1000, independent of redshift. The radiative lifetime of such electrons at 
z <~ 6 is ~ 1.0 x 10 7 years. For comparison, observations of a z — 6 jet at 1.4 GHz probe electrons with 
7 e ~ 7000, assuming a 50fiG magnetic field, corresponding to a radiative lifetime ~ 1.3 x 10 6 years. If 
particle acceleration ceased between 10 6 and 10 7 years ago, then the exponential cut-off at high energies in 
the relativistic particle population might lead to the situation observed, ic. an X-ray-loud but radio-quiet 

2 Shioya et al. (2003) have found a faint optical galaxy located within l'/9 of the quasar. They suggest that 
the QSO may be gravitationally magnified by this galaxy by a factor of about two. We have not corrected 
for this magnification in the calculations above. 
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jet. More sensitive X-ray observations, and lower frequency radio observations, are required to test this 
interesting possibility. 



5.4. Over-densities of radio sources toward the highest redshift QSOs 

For the highest redshift QSO in our sample, J1030+0524 at z = 6.28 (Figure 2), we find four fairly 
bright radio sources (Si. 4 > 200/iJy) within 1' of the QSO, one of which is an arcmin-scalc powerful double 
(FRII) radio source. Note that for a 2' field one expects only 0.5 sources by chance with Si. 4 > 200/iJy. It 
seems unlikely that the sources are at the redshift of the QSO, since an arcmin-scale double radio galaxy 
has never been detected beyond z ~ 3 (Carilli et al. 1997). Interestingly, the highest redshift QSO known, 
J1148+5251 at z = 6.4 (Fan et al. 2003), also has two bright radio sources (8 and 70 mJy at 1.4 GHz) 
within 1' of the QSO position (Bcrtoldi et al. 2003). 

A possible explanation for the excess radio source density in the fields of the two highest redshift QSOs 
is that the radio sources are in a group or cluster that happens to lie along the line-of-sight. If this is the 
case, then gravitational lensing by the cluster could magnify the QSO emission by a factor 2 or so without 
giving rise to arcsecond-scale distortions in the optical images of the QSOs. 

The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under a 
cooperative agreement with the National Science Foundation. M.A.S. acknowledges support of NSF grant 
AST-0071091. 
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Table 1. Properties of Observed QSOs 



Optical Position (J2000) S1.4 S 250 



oso 

V: ,/ kJ V / 




iviB 


RA (h m H\ 

A ViT- \l I III Uj) 






\m Ivl 


J0231-0729 a 


5.41 


-27.37 


02 31 37.65 


-07 28 54.5 


< 50 


< 3.5 


J0301+0020 b 


5.50 


-24.00 


03 01 17.01 


00 20 26.0 


73 ± 18 


0.87 ± 0.2 


J0756+4104 a 


5.09 


-26.50 


07 56 18.14 


41 04 08.6 


65 ± 17 


5.5±0.5 


J0836+0054 c 


5.82 


-28.10 


08 36 43.85 


00 54 53.3 


1750 ± 40 


< 2.9 


J0913+5919 a 


5.11 


-26.20 


09 13 16.56 


59 19 21.5 


18950 ± 380 


< 2.8 


J1030+0524 c 


6.28 


-27.37 


10 30 27.10 


05 24 55.0 


< 61 


< 3.4 


J1044-0125 d 


5.73 


-27.63 


10 44 33.04 


-01 21 49.6 


< 79 


3.4±l.l c 


J1204-0021 f 


5.03 


-27.64 


12 04 41.73 


-00 21 49.6 


< 87 


NA 


J1208+0010 s 


5.27 


-26.30 


12 08 23.82 


00 10 27.7 


< 60 


< 3.1 


J1306+0356 c 


5.99 


-27.41 


13 06 08.26 


03 56 26.3 


< 53 


< 3.1 



a Anderson et al. 2001 
b Stern ct al. 2000b 
c Fan et al. 2001 

d Djorgovski et al. 2001; Fan et al. 2000b 
c Extrapolated from 350 GHz (section 3) 

f Fan ct al. 2000a 
g Zheng et al. 2000 
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Figure Captions 

Figure 1: Images at 1.4 GHz of the ten high redshift quasars discussed in this paper. The FWHM 
of the Gaussian restoring beams are shown in the insets in all frames. Contour levels (solid lines) are a 
geometric progression in the square root of two starting at 2a, with a listed below (er corresponds to the 
measured rms on the image). Two negative contours (dashed) are included. The central cross in each image 
marks the optical QSO location. 

Figla = J0231-0728, a = 16.6 ^Jy beam" 1 ; Figlb = RD J0301+0020, a = 17.2 ^Jy beam" 1 ; Figlc = 
J0756+4104, a = 17 ^Jy beam" 1 ; Figld = J0836+0054, a = 21.6 ^Jy beam" 1 ; Figlc= J0913+5919, a = 
14 ^Jy beam" 1 ; Figlf= J1030+0524, a = 20.3 fi3y beam" 1 ; Figlg= J1044-0125, a = 26.5 ^Jy beam" 1 ; 
Figlh= J1204-0021, a = 29 ^Jy beam -1 ; Figli= J1208+0010, a = 20 ^Jy beam" 1 ; Figlj= J1306+0356, a 
— 17.8 /iJy beam -1 ; 

Figure 2: A wider field image of J1030+0524 at 1.4 GHz. The FWHM of the Gaussian restoring 
beam is 3" x 3" . Contour levels are a geometric progression in the square root two starting at 0.13 mJy 
beam -1 . Three negative contours (dashed) are included. The cross shows the position of the optical QSO. 

Figure 3: The relationship between redshift and the observed spectral index between 250 and 1.4 GHz 
for star- forming galaxies (solid curve), as derived from the models presented in Carilli & Yun (2000). The 
dashed lines show the rms scatter in the distribution. The 250-to-1.4 GHz spectral indexes for QSOs 
detected at 250 GHz are shown with la error bars. 
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Fig. 1.- 



Uppcr left: figure la, upper right: figure lb, lower left: figure lc, lower right: figure Id 
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Fig. lcont.. — Upper left: figure le, upper right: figure If, lower left figure lg, lower right figure lh 



- 23 - 




Fig. lcont.. — Left: figure li, Right: figure lj 
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